A versatile modelling approach to determine the hydrophobicity of peptides at the atomic level

This study describes a versatile computational method to determine the hydrophobicity of small peptides at the atomic level. Free energies of transfer for individual atoms in peptide structures were derived, utilising two specifically defined parameters: (i) the water-excluding distance to define the dynamic interface between a peptide solute and its surrounding solvent and (ii) the corresponding hydrophobicity index as a relative measure for water occlusion/repulsion. The method was tested on a range of small peptide models (Ac-X-NH₂, G-X-G, Ac-WL-X-LL and Ac-GG-X-GG-NH₂) and several derivatives of these structures, whereby X was any of the 20 most common amino acids that naturally occur in polypeptides or proteins. The advantage of this new method lies in its versatility, ease to implement and capability to provide information on the hydrophobicity characteristics at the atomic level. The approach also encapsulates the impact of factors that influence these properties, but which have hitherto been difficult to accurately quantify, e.g. steric hindrance or proximity effects due to nearby polarised atoms. The method is not conditional on the knowledge of hydrophobicity parameters from the literature and does not require a sophisticated computer software/hardware to enable the atomic solvent-accessible surface areas or other hydrophobicity parameters to be de novo obtained.     

The systematic approach to describing conformational rearrangements in G-quadruplexes

Conformational changes in DNA G-quadruplex (GQ)-forming regions affect genome function and, thus, compose an interesting research topic. Computer modelling may yield insight into quadruplex folding and rearrangement, particularly molecular dynamics simulations. Here, we show that specific parameters, which are distinct from those commonly used in DNA conformational analyses, must be introduced for adequate interpretation and, most importantly, convenient visual representation of the quadruplex modelling results. We report a set of parameters that comprehensively and systematically describe GQ geometry in dynamics. The parameters include those related to quartet planarity, quadruplex twist, and quartet stacking; they are used to quantitatively characterise various types of quadruplexes and rearrangements, such as quartet distortion/disruption or deviation/bulging of a single nucleotide from the quartet plane. Our approach to describing conformational changes in quadruplexes using the new parameters is exemplified by telomeric quadruplex rearrangement, and the benefits of applying this approach to analyse other structures are discussed.     

AFLP characterization of natural populations of Berberis (Berberidaceae) in Patagonia, Argentina

 AFLP markers were used to analyse the intra- and interspecific relationships among 22 natural populations of 13 Patagonian species of Berberis and the relationships among the taxa belonging to homoploid and polyploid complexes. Seven primer combinations gave rise to 231 AFLP bands, of which 199 were polymorphic. Correspondence between AFLP data, morphological traits and seed protein bands was also assessed. The dendrogram inferred from AFLP fingerprints showed that, in general, populations of the same species formed closely related groups with high coefficients of similarity. Principal co-ordinates analysis showed two separate subgroups: (i) B. bidentata and their putative ancestors –B. darwinii and B. linearifolia– which form a homogamic group, and (ii) B. buxifolia, B. heterophylla and B. parodii– which could form a polyploid hybrid complex. Received March 21, 2001 Accepted September 11, 2001     

Structural investigations into the binding mode of novel neolignans Cmp10 and Cmp19 microtubule stabilizers by in silico molecular docking,molecular dynamics,and binding free energy calculations

Microtubule stabilizers provide an important mode of treatment via mitotic cell arrest of cancer cells. Recently, we reported two novel neolignans derivatives Cmp10 and Cmp19 showing anticancer activity and working as microtubule stabilizers at micromolar concentrations. In this study, we have explored the binding site, mode of binding, and stabilization by two novel microtubule stabilizers Cmp10 and Cmp19 using in silico molecular docking, molecular dynamics (MD) simulation, and binding free energy calculations. Molecular docking studies were performed to explore the β-tubulin binding site of Cmp10 and Cmp19. Further, MD simulations were used to probe the β-tubulin stabilization mechanism by Cmp10 and Cmp19. Binding affinity was also compared for Cmp10 and Cmp19 using binding free energy calculations. Our docking results revealed that both the compounds bind at Ptxl binding site in β-tubulin. MD simulation studies showed that Cmp10 and Cmp19 binding stabilizes M-loop (Phe272-Val288) residues of β-tubulin and prevent its dynamics, leading to a better packing between α and β subunits from adjacent tubulin dimers. In addition, His229, Ser280 and Gln281, and Arg278, Thr276, and Ser232 were found to be the key amino acid residues forming H-bonds with Cmp10 and Cmp19, respectively. Consequently, binding free energy calculations indicated that Cmp10 (?113.655 kJ/mol) had better binding compared to Cmp19 (?95.216 kJ/mol). This study provides useful insight for better understanding of the binding mechanism of Cmp10 and Cmp19 and will be helpful in designing novel microtubule stabilizers.     

Molecular modelling of ionic liquids in the ordered mesoporous carbon CMK-5

We performed classical molecular dynamics simulations of the ionic liquids (ILs) [dmim⁺][Cl^?] and [emim⁺][NTf₂^?], confined in a model CMK-5 material, which consists of amorphous carbon nanopipes (ACNPs) arranged in a hexagonal array. We compare our findings against the behaviour of the same ILs inside an isolated ACNP (i.e. no IL adsorbed on the outer surface of the ACNP) and inside a model CMK-3 material (which is similar to CMK-5, but is formed by amorphous carbon nanorods). Our results indicate that the presence of IL adsorbed in the outer surface of an uncharged ACNP in CMK-5 affects the dynamics and the density of an IL adsorbed inside the ACNP and vice versa. ILs adsorbed outside the nanopipes in CMK-5 (i.e. with IL also adsorbed inside the nanopipes) have faster dynamics and remain closer to the carbon surfaces when compared to the same ILs adsorbed on CMK-3 materials. The trends are IL-specific: [dmim⁺][Cl^?] has slower dynamics when inside an isolated ACNP than when inside the ACNPs in CMK-5, but in contrast, [emim⁺][NTf₂^?] moves faster when it is inside an isolated ACNP than when it is inside the ACNPs in CMK-5 (i.e. with IL adsorbed outside the nanopipes).     

Identification of abelson tyrosine kinase inhibitors as potential therapeutics for Alzheimer’s disease using multiple e-pharmacophore modeling and molecular dynamics

Efforts to combat Alzheimer’s disease are focused predominantly on inhibiting the activity of the enzyme(s) that have been identified to be responsible for the production of the amyloid-forming peptide. However, the inherent complexity associated with the network of pathways leading to the disease may involve additional targets for designing effective therapies. Recent experimental findings have identified abelson tyrosine kinase, a non-receptor kinase as a new target for Alzheimer’s. In this work, we employed energy optimized multiple pharmacophore modeling strategy from multiple c-Abl structures bound with ligands in the inactive ATP binding conformation (DFG-out). Virtual screening followed by docking of molecules from ChemBridge resulted in the identification of 10 best scoring molecules. MD simulations of the top three complexes revealed that Compound A, C are the most stable complexes with the most persistent protein–ligand interactions consistent with the calculated binding affinities for the top three compounds. Given the implied role of c-Abl not only in AD but in Parkinson’s disease, the identified compounds may serve as leads for effective neurotherapeutics.     

Adsorption and immobilisation of human insulin on graphene monoxide,silicon carbide and boron nitride nanosheets investigated by molecular dynamics simulation

The adsorption and immobilisation of human insulin onto the bio-compatible nanosheets including graphene monoxide, silicon carbide and boron nitride nanosheets were studied by molecular dynamics simulation at the temperature of 310 K. After equilibration, heating and 100 ns production molecular dynamic runs, it was found that the insulin was adsorbed and immobilised onto the considered surfaces in a native folded state. The structural parameters, including root-mean-square deviation and fluctuation, surface accessible solvent area, radius of gyration (R_g) and the distance between the centre of the mass of immobilised protein and the surface of the considered nanosheets, were measured, analysed and discussed. The energetics of the studied systems such as the interaction energy between protein and nanosheet was also measured and addressed. The discussions were centred on the structural and energetic parameters of the protein and nanosheets, including charge density, hydrophobicity, hydrophilicity and residue polarity. The results also showed that the active site of C-termini of chain B played an important role in the adsorption process and this could be helpful in the protection of insulin in its smart delivery and release applications.     

Biomedical Enhancements as Justice

Biomedical enhancements, the applications of medical technology to make better those who are neither ill nor deficient, have made great strides in the past few decades. Using Amartya Sen's capability approach as my framework, I argue in this article that far from being simply permissible, we have a prima facie moral obligation to use these new developments for the end goal of promoting social justice. In terms of both range and magnitude, the use of biomedical enhancements will mark a radical advance in how we compensate the most disadvantaged members of society.